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I •  Project  Summary 

Ibis  report  summarizes  the  work  undertaken  between  1  July 
1973  and  1  January  1974,  and  should  be  considered  as  a  supple¬ 
ment  to  the  previous  report. 


II •  Technical  Report 


The  principal  effort  for  this  period  has  been  the  develop¬ 
ment  of  computer  programs  for  the  processing  of  high  density 
data  acquired  at  sea  and  preparation  for  a  cruise  to  the  central 
North  Pacific  in  February-March  1974. 

A  program  which  deconvolves  the  raw  digital  data  record 
from  the  two  temperature  gradient  channels  and  the  conductivity 
channel  has  been  debugged  and  is  now  in  service.  This  permits 
combining  the  low  and  high  frequency  temperature  (or  conductiv¬ 
ity)  information  from  two  separate  data  channels  into  one  record 
of  temperature  vs.  depth.  The  dynamic  range  achievable  with  this 
technique  is  greater  than  105,  thus  permitting  a  single  record  to 
contain  virtually  all  of  the  significant  information  on  the  ver¬ 
tical  profile.  The  resultant  data  is  being  analyzed  to  determine 
the  patch  sizes  of  regions  of  intense  microstructure  and  their 
relation  to  large-scale  irregularities. 

During  the  period  21  February  -  20  March  1974,  vertical 
profiles  will  be  ma  e  at  28°N,  155°W  in  the  central  North  Pacific, 
making  the  third  set  of  data  from  this  location.  The  other 
measurements,  made  in  September-October  1971,  and  June  1973,  will 
be  compared  with  this  set  in  an  attempt  to  determine  seasonal 
variations  in  levels  of  microstructure  activity  in  the  upper  two 
kilometers  and  to  relate  these  to  changes  in  the  heat  content  of 
the  wate~  column. 
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I. 


Project  Summary 


This  report  supplements  previous  reports  on  development  and 
use  of  a  freely  floating  instrumented  capsule.  A  summary  of 
work  undertaken  July  1,  1973  throi  ’’h  December  31,  1973  is  given. 
Future  plans  are  indicated. 


II .  Technical  Report 

A.  Our  major  effort  during  this  period  has  been  on  the  analysis 
of  data  gathered  in  June,  1973.  At  that  time  we  participated  in 
a  joint  OWEX  experiment  with  other  investigators  in  an  area  300 
miles  S.W.  of  San  Diego.  We  obtained  three  thermal  microstruc¬ 
ture  records,  each  of  one  day  duration,  and  one  internal  wave 
record  of  6  days  length. 

1.  Analysis  of  the  intei nal  wave  data  is  essentially  com¬ 
pleted  and  will  be  reported  during  the  early  part  of  1974. 

Figure  1  gives  a  preliminary  look  at  two  internal  wave  spectra 
each  estimated  from  half  of  the  six  day  internal  wave  record. 

The  observed  spectra  decrease  as  a;-2  up  to  the  local  bulk 
vaisala  frequency  and  drop  off  sharply  at  higher  frequencies. 

Detailed  thermal  structure  profiles  typically  reveal  thermal 
inversions.  In  the  present  experiment  we  obtain  such  profiles 
repeatedly  while  advecting  with  the  mean  water  motion.  By  remov¬ 
ing  the  internal  wave  oscillations  from  the  records,  it  is  possi¬ 
ble  to  more  clearly  observe  changes  in  the  thermal  structure. 
Figure  2  is  a  preliminary  sketch  of  a  five  hour  segment  of  a 
structure  record  referenced  somewhat  imperfectly  to  the  4.665°C 
isotherm.  It  shows  a  thermal  inversion  pattern  suggestive  of  a 
Kelvin-Helmholtz  shear  instability.  Four  inversions  of  this  type, 
each  lasting  several  hours,  were  identified  in  a  three  day  segment 
of  the  data. 

2.  Several  inversions  of  similar  scale  appear  in  the  records 
of  drops  which  included  instrumentation  for  measuring  thermal 
microstructure.  The  microstructure  "signatures”  of  these  inver¬ 
sions  will  be  used  to  help  understand  their  physics.  Kelvin- 
Helmholtz  shear  instability,  already  mentioned,  or  interfingering 
of  statically  stable  water  masses  with  salinity  compensated  tem¬ 
perature  inversions  are  two  possibilities.  The  analysis  of  all  of 
the  June  data  is  expected  to  be  completed  by  mid  1974. 

B.  A  two-axis  self-contained  electromagnetic  current  meter 
designed  and  prepared  by  Jack  Olson  of  NUC  was  included  in  the 
midwater  float  instrumentation  m  initial  sea  trials  in  November. 
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III.  Future  Plans 


A.  Modifications  to  the  midwater  float  and  portable  laboratory 
are  underway  in  preparation  for  a  nominal  30  day  internal  wave 
measurement  to  be  conducted  mid- June  to  mid-July. 

B.  Improvements  in  the  electromagnetic  current  meter  based  on 
results  of  the  November  sea  trials  are  being  made  and  a  second 
unit  is  being  built  for  shear  measurements  in  conjunction  with 
microstructure  and  internal  wave  measurements. 

C.  Coordination  and  planning  with  Mike  Gregg  has  begun  in  pre¬ 
paration  for  joint  microstructure  and  internal  wave  measurements 
in  the  fall. 
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The  solid  curve  is  the  spectrum  estimated  from  the  first 

dashed dayS  of.the,  S1X  daV  isotherm  depth  series.  The 

A  s£rfiqhtVli^  £5®  ?Stima-f  fr°m  the  second  th^e  days. 

^  ^ne  of  sl°Pe  "  is  given  for  reference.  N 
denotes  the  local  Brunt-Vaisaia  frequency. 
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Figure  2. 


Ihis  figure  shows  a  five  hour  segment  of  the  thermal 
i^rS?Sre  motl<?n  relative  to  an  isotherm.  The  ordinate 
reJatlVe  to  that  isotherm.  The  heavy  curve 
are  contours  of  isotherms  spaced  every  .005°C,  and  the 
various  shaded  areas  are  areas  where  the  temperature 
lies  within  particular  .005°C  bands,  as  indicated  by  th 

takljfbv  the  Jight2r  sawtooth  curve  is  the  path  in  time 
taKen  by  the  temperature  sensor. 
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Project  Summa^ 


As  described  in  the  previous  Technical  Report,  a  well 
instrumented  open  ocean  experiment  OWAX  was  successfully  under- 

ofkthe  r^ln9  t^e/irS,t  Six  months  of  19?3.  The  primaryYgoals 
31  1973  f°r  the  Period  from  July  1,  1973  to  December 

inertial  «.  V  t0  derive  equations  to  utilize  the 

dSS  ?ata  to  correct  the  velocity  and  wave  height 

of^V  motion ,  2)  to  undertake  preliminary  analysis 

of  the  atmospheric-turbulent  boundary  layer  data,  and  3)  to 

data^o  atmos!?^ric  boundary  layer,  internal  and  surface  wave 
is  proceeding  correla tions ■  Detailed  analysis  of  the  data 


11 •  Preliminary  Results 

--atform  Correction  for  FLIP'S  Motion.  Measure- 
tuJbnipnf  fe  vaFious  velocity  components  in  the  atmospheric 
fc  ^ndary  layer  as  well  as  wave  height  information 

correotJ^fF^n  3  non-lnertial  platform  such  as  FLIP  must  be 
orrected  for  the  motion  of  the  platform  or  sensors.  This 

-irs^thr^T1^^?  Ue  measured  data  in  two  "lain  respects: 

level  vatiot  h®  Pr°beS  ab°VG  the  mean  sea  surface 

citv  ic k ltJ,tlme'  and  thus  a  spurious  time  dependent  velo- 
*  ead  bY  the  anemoineters  because  of  the  presence  of  the 
velocity  Profile  in  the  atmospheric  boundary  layer 

an  Unknown  !h°  themselves  are  being  swept  through  space  with 
3*  !!  fc  varying  velocity  with  respect  to  an  inertial 

r lystem  and  this  velocity  must  be  subtracted  from  the 
recorded  data.  During  OWAX,  FLIP's  motions  were  accurately 

soCthatdihY  mGanS  °f  3  Litton  inertial  platform  and  gyroscopes 
^hat.fche  anemometer  data  could  be  corrected.  The  f  lowing 

j  briff  description  of  a  correction  procedure  we  developed 
and  which  will  be  applied  to  the  data  in  the  near-future. 

g3ne3aiJ motion  of  a  rigid  body  (FLIP)  can  be  considered 
as  composed  of  the  translational  motion  of  any  point  in  that  body 

no  a  rigid  body  rotation  about  that  point.  On  FLIP  the  trans- 

«ad  Sv  °btaine?  b?  int^atiSn  of  the  acceferaUoS 

ed  from  t-hp  ^htt0n  ln?rtlal  Platform,  and  the  rotation  is  obtain- 
from  the  three  angles  read  from  the  gyroscopes.  Since  the 

gyroscope  axes  maintain  their  orientatio^  in  space,  this  oriSnta- 

comSn1^  ?Sed  the  direction  of  the  desired  velocity 

opponents,  with  the  z  axis  (yaw)  being  vertical  to  the  mean 

wind11  3nd  thS  X  axis  Pointing  approximately  into  the 

nd.  Let  the  gyroscope  coordinate  system  be  denoted  by  S,  and 

the  system  fixed  to  FLIP  be  denoted  by  s.  Then  it  can  be  shown 
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that. any  point  fixed  to  FLIP  /described 
as  r1)  will  have  components  r1  in  S 
given  by  the  following  transformation: 
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Where  the  angles  in  the  rotation  matrix  are  with  respect  to  the 
gyroscope  gimbaling  system  and  are  read  out  directly  as  data: 


0  =  yaw  (z  axis) 

=  pitch  {approx,  y  axis) 
ip  =  roll  (approx,  x  axis) 


n  particular  the  fixed  values  r  are  chosen  to  be  the  vector 
components  of  the  probe  location  with  respect  to  the  Litton 
inertial  platform  at  time  t  =  0  when  e  =  <j>  =  y  =  o,  and  these 
are  known  from  the  geometrical  placement  of  the  apparatus. 

The  total  motion  of  the  probe  with  respect  to  an  inertial 
coordinate  system  can  then  be  given  in  two  parts: 
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vj  (0) 
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+  /  a](t')dt' 


integration  of  iner¬ 
tial  platform  aata 


{—  (t)  }rk 

dt  k 

^ _  ,  - 


derivative  of 
gyroscope  data 


(2) 


j  =  1/  2,  3 

The  probe  position  as  a  function  of  time  can  be  obtained 
by  further  integration: 

t  t' 

x3(t)  =  (0)  +  v3(0)t  +  J  J  a  3  (t " )  dt "  + 

p  a  a  0  0 


(t)ry 
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(3) 


The  two  constants  of  integration  are  not  explicitly  given 
by  the  experiment  but  mav  be  approximately  found  by  the  follow¬ 
ing  argument.  If  FLIP  is  adequately  moored,  one  expects  its 
time  average  motion  to  vanish.  Hence  by  letting 

(t)  =  0  and  x^  (t)  =  0 , 

P  p 

one  obtains  the  initial  values  v^(0)  and  x^(0). 

a.  a 

By  the  above  technique  the  instantaneous  location  and 
velocity  of  the  probe  is  known,  and  thereby  the  instantaneous 
velocity  data  can  be  corrected: 
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It  must  be  further  noted  that  since  the  velocity  data  is  read 
by  sensors  which  change  orientation  with  FLIP'S  motion,  the 
interpretation  of  the  velocity  components  needs  to  be  corrected 
by  the  same  rotation  matrix  that  was  exhibited  above.  Then 
finally: 
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In  light  of  the  above  discussion  several  errors  arising  from 
FLIP's  motion  need  to  be  determined  numerically  by  computer  cal-- 
culation.  By  use  of  equation  (3)  the  vertical  location  of  the 
probe  is  computed  as  a  function  of  time.  Assuming  that  the  meas¬ 
urements  are  made  in  the  logarithmic  portion  of  the  atmospheric 
boundary  layer,  it  is  possible  to  estimate  the  error  in  the  meas¬ 
ured  signal  due  to  the  changing  vertical  location  of  the  probe. 

The  three  correction  terms  appearing  in  equation  (5)  need  to  be 
computed  and  separately  assessed.  By  cross-multiplication  of 
equation  (5)  and  time  averaging,  one  may  directly  assess  the 
relation  between  the  corrected  and  uncorrected  Reynolds  stresses. 
By  taking  the  Fourier  transform  of  each  term  separately,  one  may 
see  if  the  correction  terms  add  a  dominant  spurious  spike  into 
the  spectra  of  fluid  turbulence.  Another  effect  of  probe  motion 
which  needs  to  be  considered,  but  which  one  expects  to  be  negli¬ 
gible,  is  the  apparent  broadening  of  the  frequency-wavenumber 
relation  that  is  customarily  obtained  from  Taylor's  hypothesis. 

In  other  words,  a  probe  moving  into  the  wind  records  higher  fre¬ 
quency  turbulence  than  while  moving  away  from  it.  This  effect  is 
felt  in  scalar  spectra  as  well  as  velocity  spectra. 

B*  Atmospheric  Turbulent  Boundary  Layer  Data.  Analysis  of  OWAX 
temperature  and  humidity  data  has  confirmed  other  evidence  indi¬ 
cating  that  the  temperature  spectrum  shows  a  definite  departure 
from  the  -5/3  law  inertial  subrange  form  expected  for  passive 
scalar  mixing  at  high  Reynolds  number.  The  humidity  spectra 
exhibit  the  expected  inertial  subrange  behavior  as  shown  in  Fig¬ 
ure  11  of  the  previous  report.  It  is  interesting  to  note  that 
using  the  same  electronic  equipment  and  techniques,  measurements 
of  temperature  and  humidity  spectra  in  the  atmospheric  boundary 
layer  over  land  were  carried  out  this  summer  at  a  Minnesota  test 
site  in  a  joint  experiment  with  the  Air  Force  Cambridge  Research 
Laboratories.  Both  temperature  and  humidity  spectra  exhibited 
the  expected  -5/3  behavior  over  land.  The  reasons  for  this 
result  is  still  under  investigation  and  seem  of  fundamental  impor¬ 
tance  to  the  physics  of  the  marine  boundary  layer. 

C .  Internal  Wave,  Surface  Wave  and  Atmospheric  Boundary  Layer 
Interactions .  Some  preliminary  indication  exists  that  internal 
waves  are  influencing  sea  surface  temperature  and  atmospheric 
temperature.  From  some  strip  chart  data  recordings,  the  air 
temperature  can  be  seen  to  decrease  and  the  sea  surface  tempera¬ 
ture  increases  as  an  internal  wave  event  occurs.  The  cause  and 
magnitude  of  the  surface  response  to  the  internal  wave  event,  as 
well  as  exactly  what  is  measured  by  the  Pinkel  detector,  are 
under  study. 
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